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Level spacings of parametric chiral random matrices 
and two-color QCD with twisted boundary condition 

Shinsuke M. Nishigaki 
Graduate School of Science & Engineering, Shimane University, Matsue 690, Japan 

We evaluate level spacing and smallest eigenvalue distributions of chiral random ma- 
trix ensembles transiting from symplectic or orthogonal to unitary symmetry classes with 
a crossover parameter p. As expected from the effective a model description, these results 
can be fitted perfectly to the fundamental or adjoint staggered Dirac spectrum of SU(2) 
quenched lattice gauge theory under the imaginary chemical potential (twisting) fj,. The 
linear dependence of the parameter p on p determines the pion decay constant F as its 
proportionality constant. 



Subject Index: 100, 138, 160, 169, 183, 350 

Introduction. Two-color QCD has served as a strategic testing ground for the realistic 
chromodynamics, as well as a tractable model interesting by its own right, due to the 
absence of the sign problem at finite density and with pairs of degenerated flavors0 
This and other gauge theories with quarks in (pseudo)real representation exhibit an 
exotic global symmetry breaking® quarks and charge-conjugated antiquarks form a 
i— i multiplet of the extended flavor group, which breaks spontaneously to the extended 

vector subgroup© Accordingly the effect of the chemical potential that distinguishes 
quarks from antiquarks and breaks this extended symmetry is incorporated into 
the low-energy chiral Lagrangiar® through the flavor-covariant derivative BE As 
the unconventional global symmetries originate from the presence of antiunitary 
symmetries of Dirac operators, the symmetry- violating chemical potential in these 
p<| QCD-like theories manifests itself also in the statistical properties of Dirac spectral^ 1 

One can indeed predict the fluctuation of Dirac eigenvalues that permeate into the 
complex plane from the zero- momentum part of chiral Lagrangians, which in turn is 
equivalent to chiral Gaussian Orthogonal or Symplectic Ensemble (chGOE, chGSE) 8 ' 
extended to its non-Hermitian counterpart through the introduction of a schematic 
OO density component coupled to (real) /xEMlH} 

Antiunitary symmetries of SU(2) Dirac operators could, however, as well be 
violated by any Hermitian term in complex representations, most simply by the 
U(l) gauge field either fluctuating or fixed as a background, the latter beinggauge- 
equivalent to the twisted boundary condition or imaginary chemical potentialllll'llHf'lHf 
In this case Dirac eigenvalues stay real after the inclusion of symmetry violation and 
their statistical behavior exhibits crossover, instead of permeation to the complex 
plane. Since these two cases are essentially identical in the chiral Lagrangian descrip- 
tion save for the sign of [i 2 , studies of Dirac spectra from both sides would reveal 
the validity of analytic continuation in /j>, as much required for three-color QCD. 

Spectral crossover between universality classes of Hermitian random matrix 
(RM) ensemblesJU'Ei'OIJ namely from GOE or GSE to GUE, has been extensively 
studied for disordered and chaotic Hamiltonians whose time-reversal invariance is 
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slightly broken by a magnetic field. The transitional behavior of spectral fluctua- 
tions is by itself universal, in a sense that the local spectral fluctuation depends only 
on a single parameter p defined below. The reason for this universality is traced 
back to the nonlinear a model (eq.Q below) governing the spectral statistics, which 
can either be derived by conventional disorder averaginglil 1 or by summation over 
encountering periodic orbits of chaotic systemsJlH' without any reference to the de- 
tails of dynamics. Thus one can rely on the simplest model that yields the identical 
a model, i.e. parametric RMs, for actual computation. 

On the other hand, chiral or superconducting variants of universality crossover 
are relatively less explored. Damgaard and collaborator£U''[l2J have marked a signif- 
icant breakthrough by analytically computing correlation functions and individual 
small eigenvalue distributions for the spectral crossover within the chiral Gaussian 
Unitary Ensemble (chGUE) class. They presented a convincing numerical evidence 
that the Dirac spectrum of three-color QCD at imaginary isospin chemical potential 
indeed exhibits crossover as predicted by chiral RMs. On the other hand, although 
analytic results for microscopic spectral correlation functions for the crossover from 
chGOE or chGSE to chGUE are known for some time Jill 1 JHI they have lacked phys- 
ical application. To the best of our knowledge, the only physical example of the 
crossover involving different Hermitian chiral universality classes is the CI-C transi- 
tion for the super /normal/super-conducting hybrid interface in a magnetic field© 
In this Letter we provide novel applications of crossover between chiral Hermitian 
universality classes from lattice gauge theory. 



Parametric chiral random matrices. Consider an ensemble of iV x N Hermitian 
complex (quaternion) matrices H = Hs+iaH^, with H$ real symmetric (quaternion 
selfdual) and H\ real antisymmetric (quaternion antiselfdual) , distributed according 
to Gaussian measures of variance a 2 . By setting N/2x N/2 block-diagonal parts of H 
to zero, the matrix takes the form ( T and D stand for transpose and quaternion-dual) 



H 



Hi+ iaH 2 

(Hi - iaH 2 ) T > D 



N N 

Hi, H 2 : — x — (quaternion-)real matrix. 



Depending on the parameter a, this parametric (also called dynamical or Brownian- 
motiorllll') chiral RM ensemble interpolates between two limiting cases, chGOE 
(chGSE) at a = and chGUE at a = 1. Since nonzero eigenvalues of H occur 
in pairs of equal magnitude and opposite signs, if suffices to retain non- negative 
eigenvalues only. The correlation function of n eigenvalues {Aj} of H in the vicin- 
ity of the origin (where the mean level spacing is A(0)) is expressed, in the limit 
N —7- oo, a —7- and p = ao~/A(Q) fixed, as a determinant with S, D, I given by^ ; ' ~~ 1 

S(x,y) I(x,y) 
D(x,y) S(y,x) 



R n (xi, . . . ,x n ) = \/det[K(xi,Xj)]" j=1 , K(x,y) 



S{x, y) = iTy/xy { ^WW-^oW^N + M^y) j ^ e - P H^)j o{vx) 

x y 2 



D(x,y) = — / dvv 2 e 2pv {xJi(vx)J Q (vy) - J (vx)yJi(vy)} , 
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I(x,y) = ^ J°° dvv J°° due'^ 1 ^ {J (vux)J (vy) - J (vx)J (vuy)} (2) 
for chGOE-chGUE crossover, and 



x 2 — y 2 



s(x,y) = n ^l *M™)JoW-M™)vJiW _ M^l f dv ^i^) Mvy) ) f 

Jxy 

D{x,y) = 

I ( x ,y) = ^Y- J dvv 2 e~ 2pv {xJi(vx)J (vy) - J (vx)yJ 1 (vy)} (3) 

for chGSE-chGUE. Here Xi = Aj/Z\(0) are unfolded eigenvalues. The probability 
E(s) that the interval [0, s] contains no eigenvalue is given as a Fredholm determinant 







dv v 


/ due' 




'0 


dv v 2 


e~ 2 ^ 2 



E(s) = •\/Det(l - K s ), (4) 

where K s is an integral operator of convoluting with the dynamical Bessel kernel 
K(x,y) ([l])~(p]) acting on two-component L 2 -functions / over the interval [0, a], 
(K s f)(x) = f Q dy K(x,y)f(y). Probability distribution p\(s) of the unfolded small- 
est eigenvalue s = Xi/A(0), which is a half of the very central level spacing, is given 
by the first derivative pi(s) = —E'(s). SedHf for an alternative derivation. 

In the limit i,i/>1 with x — y kept finite, eigenvalues x, y are liberated from 
repulsion by their mirror images — x, —y and their correlation becomes translationally 
invariant. The S,D,I components of the kernel in this limit are given by, as functions 
of r = x — y^i 

S(r) = , D(r) = / — ve 2p v sinw, J(r) = / — e~ 2p v sinwr, (5) 

irr In ir .L ttv 



sin7rr f°° dv _ 2o 2 v 2 ■ n ^ H dv 2 ^ a . 

S(r) = , D(r) = / — ve p smvr, I(r) = / — e p smvr, (6) 

vrr J n vr J TTV 



for GOE-GUE and GSE-GUE crossover, respectively. These expressions interpolate 
between two nonchiral ensembles, GOE (GSE) at p = and GUE at p = oo. The 
gap probability E(s) is again given as Q with this dynamical sine kernelEU 1 The 
distribution P(s) of level spacings s = Xi + \ — X{ is given by the second derivative 
P(s) = E"(s). 

An efficient way of evaluating Fredholm determinant of a trace-class operator K s 
acting on L 2 -functions over an interval [0, s] is the Nystrom-type discretization!!} 'Hii 

Det(l - K s ) ~ det [Jy - Kfa, Xj)y/WiWj]™ j=1 . (7) 

Here the quadrature rule consists of a set of points {x%} taken from the interval [0, s] 
and associated weights {w;} such that J Q S f(x)dx ~ X^2=i f( x i) w i- As the order m 
of the approximation increases, the RHS of is proven to converge uniformly to 
its LHS. The convergence is rapid and exponentially- fast© In our purpose we em- 
ploy Gauss quadrature rule (sampling at the Legendre nodes), as 15-digit accuracy 
is known to be attainable already with only m = 5 for the Fredholm determinant 
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Fig. 1. Smallest eigenvalue distributions (top), microscopic level densities (center), and level spac- 
ing distributions (bottom) for (ch)GSE-(ch)GUE (left) and (ch)GOE-(ch)GUE (right) crossover. 



E(0.1) for the sine kernelH) We have applied the Nystrom-type method to the 
dynamical Bessel ([I]), Q, ^ and sine kernel ([l]), ([5]), Q, and evaluated pi(s) and 
P(s) for chG(0,S)E-chGUE and G(0,S)E-GUE crossover, respectively. In order to 
achieve accuracy needed for computing the first or second derivatives (pi(s) or P(s)) 
to a good precision, we chose the approximation order m to be at least 20 for the 
former and 100 for the latter, and confirmed the stability of the results under the 
increment of m. Plots of pi(s) and P(s) for the region < s < 3 ~ 4 and for the 
parameter range p < 1 are exhibited in Fig. 1. Although Mehta-PandejE!' have 
expressed P(s) in terms of eigenvalues of an infinite-dimensional matrix, elements of 
which are integrals involving prolate spheroidal functions, numerical plots of pi(s) 
for parametric chiral RM ensembles have never appeared in the literature. We also 
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exhibit single- level densities Ri(x) ([!]) in Fig. 1, whose first peaks are comprised of 
Pi(s). Practical advantage of adopting distributions of individual level spacings over 
n-level correlation functions (Ri(x), etc.) for fitting is clear from the figures. As the 
oscillation of the latter consists of overlapping of multiple peaks, the characteristic 
shape of each peak is inevitably smeared, only to yield a rather structureless curve 
for which an accurate fit is difficult. On the other hand, the shape of the former is 
clearly distinguishable and is extremely sensitive to the p parameter, for the ratio 
of pi(s) or P(s) for the orthogonal and symplectic classes to that for the unitary 

2 2 

class grows as exp for large s. Therefore they admit very sharp one-parameter 
fitting by the least square method or (in principle) simply from the tails of the curves. 

Dirac Spectrum. Dirac operator of QCD-like theory with quarks in real or pseu- 
doreal representation possesses antiunitary symmetry unlike those in complex rep- 
resentations HI Euclidean Dirac operator for quarks in the fundamental (adjoint) 
representation of SU(2) commutes with CT2K (CK). Here C is the charge conju- 
gation matrix, T2 is one of the generators of the gauge group, and K the complex 
conjugation. As (CT2K) 2 = +1 ((CK) 2 = —1), D is essentially a real symmetric 
(quaternion selfdual) matrix. The reality and selfduality of the continuum Dirac 
operators are known to be interchanged for the lattice staggered Dirac operators, 
due the absence of the charge conjugation matrixMf Since the Dirac operator in the 
fundamental representation of U(l) in continuum or on a lattice possesses no such 
antiunitary symmetry, Dirac operator in SU(2)-fundamental (SU(iV)-adjoint)xU(l) 
has its antiunitary symmetry (weakly) broken. As a simplest example we consider 
SU(2) quenched lattice gauge theory under the twisted boundary condition, that is 
to multiply SU(2) link variables on the temporal boundary of the hypercubic lattice 
of size V = L A by a constant phase e^"^ with n ^ = 2-Kip S^uS^ (ip *C 1). As 
SU(2) Dirac operators possess (pseudo)reality either for the periodic or antiperi- 
odic boundary condition on each dimension, we impose periodicity (antiperiodicity) 
for the spatial (temporal) direction and consider a small deviation (twisting) in the 
temporal boundary condition. This twisting is gauge-equivalent to a fixed U(l) 
background of flux 2irip or imaginary chemical potential p = i2ir<p/L, which is the 
measure of symmetry violation and plays the role of a in the parametric RMsH3 Its 
effect on the chiral Lagrangian is completely dictatedJiHf in parallel with the case of 
real chemical potentials HI 

For our aim of confirming the presence of chG(S,0)E-chGUE crossover in lattice 
gauge theories, we restrict ourselves to the strong coupling region of SU(2) at (3 = 
A/g 2 = ~ 1, where the level density at the origin, 1/Z\(0), is sufficiently above zero 
and the chiral symmetry is spontaneously broken. Accepting that the model is away 
from the continuum limit, we employ the simplest algorithm: unimproved plaquette 
action and 10-hit heat-bath update coupled with overrelaxation. Because of our need 
to detect possibly small deviations of spectral fluctuation from the universal RM 
statistics at either end (p = or 00), we give priority to the number of independent 
gauge configurations and perform our simulation on a lattice of the smallest size 
V = 4 4 . This choice is sufficient for measuring the short-separation behavior of 



6 



Letters 



eigenvalues within 3A ~ 4A and determining the p parameter precisely. In this 
region the systematic deviation due to the smallness of the lattice is expected to be 
less prominent (it shall manifest at larger separation) than the statistical fluctuation. 

Dirac spectral statistics are fitted to parametric RM predictions by the follow- 
ing steps: (i) First we perform pure SU(2) simulations for each (3 and measure the 
smallest fundamental or adjoint staggered Dirac eigenvalue Ai for O(10 5 ) configura- 
tions. Taking for granted that Dirac eigenvalues obey chG(S,0)E statistics, find the 
value of A that optimally fits the histogram of the smallest Dirac eigenvalue to the 
rescaled chG(S,0)E result p\(\\/A)/A. (ii) Next we multiply SU(2) link variables 
by the twisting phases and measure all Dirac eigenvalues {Aj} for A^onf = O(10 4 ) 
independent configurations. The unfolded smallest eigenvalue is still defined by 
%i = Ai/Z\(0) with respect to A(0) determined from pure SU(2) simulations. We fit 
the frequencies of x\ to the RM prediction pf^ (s) by the least square method and 
find the optimal value of p. The valid range of fitting the smallest eigenvalue x\ is 
chosen to be [0, 2.8] (fundamental) and [0, 1.6] (adjoint), divided into 20 segments, 
(iii) Finally we measure the distribution of level spacings from the spectral 'plateau' 
[A m , Am] adjacent to but not including the origin, in which the mean level spacing 
A{\) is well approximated as a constant close to A(0). In order to avoid possible 
distortion of the level spacing distribution, we take A m not too close to the origin and 
set it to be the 11 th eigenvalue at the smallest. We fit the frequencies of unfolded level 
spacings s = (Aj + i — Aj)/Z\(A), Aj,Aj+i G [A m ,Afvi] to the RM prediction P( p \s) as 
before. Thanks to the enormous gain in statistics by the spectral averaging,xxxx we 
can safely set the fitting range to be as large as [0, 3.8] and divide it into 40 segments. 

Simulation results. We generated 40000 independent configurations at each value of 
the coupling constant j3 = 0, 0.5, 1 and the twisting ip = 0.01 ~ .06, on a lattice of 
dimensions V = 4 4 . Optimal values of p determined from the smallest eigenvalue dis- 
tributions (SED) and level spacing distributions (LSD), for the fundamental (F) and 
adjoint (A) representation, are tabulated in the central columns of Table I. Sample 



Table I. Crossover parameters and low-energy constants 





rep/dist 


A 


E 
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\fAp 




F 2 










<^=.01 


.02 


.03 


.04 


.05 


.06 


H 


s 







F/SED 
F/LSD 


.00930(2) 
.00929(0) 


1.319(2) 


.059(1) 
.061(0) 


.121(1) 
.122(0) 


.182(1) 
.182(1) 


.239(2) 
.244(1) 


.302(2) 
.304(2) 


.363(2) 
.364(4) 


.370(1) 
.374(1) 


.215(1) 
.220(1) 


.290(1) 




A/SED 
A/LSD 


.00620(1) 
.00618(0) 


1.980(3) 


.076(1) 
.075(3) 


.150(2) 
.149(4) 


.224(3) 
.223(6) 


.288(5) 
.298(9) 


.373(7) 
.371(14) 


.435(9) 
.444(22) 


.372(3) 
.372(5) 


.217(3) 
.218(6) 


.432(11) 


0.5 


F/SED 
F/LSD 


.01014(1) 
.01004(0) 


1.210(1) 


.052(1) 
.052(0) 


.102(1) 
.103(0) 


.156(1) 
.155(1) 


.206(2) 
.206(1) 


.257(1) 
.257(1) 


.308(2) 
.307(2) 


.330(1) 
.329(1) 


.171(1) 
.170(1) 


.208(1) 




A/SED 
A/LSD 


.00629(1) 
.00622(0) 


1.951(2) 


.066(1) 
.067(3) 


.131(2) 
.134(4) 


.197(3) 
.202(5) 


.258(4) 
.268(7) 


.312(5) 
.331(9) 


.370(7) 
.400(18) 


.324(2) 
.336(4) 


.165(2) 
.177(4) 


.350(9) 


1 


F/SED 
F/LSD 


.01132(13) 
.01105(0) 


1.084(12) 


.046(2) 
.048(0) 


.094(1) 
.095(0) 


.138(1) 
.143(1) 


.185(1) 
.190(1) 


.233(2) 
.237(1) 


.280(2) 
.284(2) 


.315(1) 
.319(1) 


.156(1) 
.160(1) 


.177(1) 




A/SED 
A/LSD 


.00633(2) 
.00632(0) 


1.939(7) 


.066(2) 
.067(4) 


.134(2) 
.133(4) 


.198(3) 
.198(5) 


.257(4) 
.265(7) 


.320(6) 
.331(11) 


.391(7) 
.398(16) 


.332(2) 
.335(4) 


.173(2) 
.176(5) 


.342(9) 




Fig. 2. Smallest eigenvalue distributions (top), microscopic level densities (center), and level 
spacing distributions (bottom) of SU(2)-fundamental (left) and adjoint (right) staggered Dirac 
operators at coupling j3 = 0.5 and flux ip = 0.01 ~ .06. V = 4 4 , iV conf = 40000. 



plots of pi(s), Ri(x), and P(s) at /3 = 0.5 are exhibited in Fig. 2. In all figures, 
measured histograms are plotted by colored dots, and optimally-fitting parametric 
RM results are shown in curves of the same color. Also plotted in the figures are the 
results from chG(S,0)E or G(S,0)E (black real line), and chGUE or GUE (broken 
line). Note that the microscopic level densities Ri(x) are not fitted to the correspond- 
ing data, but are merely substituted the values of p determined from the SEDs. Even 
at inspection, the precision of one-parameter fitting is convincing in all distributions 
shown. The x 2 /d.o.f. for the optimal RM distributions are: 0.53 ~ 1.47 (Fund/SED), 
0.56 ~ 1.57 (Adj/SED), 0.64 ~ 1.49 (Fund/LSD) 0.68 ~ 1.30 (Adj/LSD). Con- 

Except for Fund/LSD at very small p < 0.1, where the distribution becomes extremely peaky 
at small s due to the onset of Kramers degeneracy and the fitting error is inevitably enhanced. 
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sidering the smallness of our lattice, the precision achieved is astonishing; it is even 
more surprising when one recalls the acute sensitivity of p( p \s) on p shown in Fig.l. 
These goodness-of-flts are comparable to those in the pioneering papers^ fitted to 
the spectral data from larger lattices: x 2 /d.o.f.= 0.33 for quenched QCD on 12 and 
X 2 /d.o.f.= 1.13 ~ 1.33 for dynamical QCD on 6 4 , with O(10 3 ) configurations. 
From Table I one notices that (i) for 
a fixed gauge coupling /?, ip-p plots are 
all quite linear. In Fig.3 we present a 
sample plot of y/Ap/ip at (5 = 1. Er- 
ror bars indicated are statistical only. 
The mean values of this ratio, namely 
from LSDs, are very stable under the 
change of flux <p and fluctuate only within 
0.5%, and (ii) the values of p deter- 
mined from SED and from LSD are in 
a good agreement, as should be. High 
linearity will be essential for the pre- 
cise determination of the pion decay con- 
stant. 




Fig. 3. Ratios between the twisting ip and 
the crossover parameter V~Ap, at /3 = 1. 



Low-Energy Constants. The effective low-energy Lagrangian for QCD-like theories 
with Np flavors of quarks in (pseudo)real representation, at finite chemical potential 
p and bare quark mass m is unambiguously fixed by the global symmetry alone (as 
long as p is much smaller than the vector meson mass) and takes the form containing 
two phenomenological constants: F the 'pion' decay constant and E = ^ip) /Np 
the chiral condensate, both measured in the chiral and zero-chemical potential limit 
m, p — > 0. If the theory is in a finite volume V = L 4 and Thouless energy defined 
as E c ~ F 2 /EL 2 is much larger than m, the path integral is dominated by the zero 
mode and takes the tractable form 

Z= [ dU exp[Vp 2 F 2 tT(BU j BU + BB) + -VIJmRetTMu) . (8) 

Js\J(2N F ) V 2 / 

Here U is an S\J(2Np) matrix, B = a-^lNp, M = ai0lN F (io~2 <X> 1jv>) f° r quarks in 
real (pseudoreal) representation. In order to extract Dirac spectrum, one introduces 
fermionic as well as bosonic quarks in the fundamental theory, leading to the graded 
group version of ^ on the effective theory side. Parametrizing the graded matrix 
U in terms of its eigenvalues and comparing the resulting expression (after analytic 
continuation p — )• ip and m — > iX) with the RM results ^ and ([3]), the coefficients 
of chemical-potential and 'mass' terms in the exponents on both sides are readily 
identified as 4VF 2 p 2 = 2n 2 p 2 and VEX = irx. The latter is merely the definition 
of the unfolded eigenvalues x = X/A because of Banks-Casher relation E = tt/AV 
that determines one of the low-energy constants E in terms of A. Eliminating the 
volume in favor of the level spacing, the former equation reads 
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where the left hand side is a volume-independent combination. Accordingly one 
can determine another low-energy constant F 2 / E from the slope of tp-yAp plots, 
preferably on lattices of various sizes. In the next-to-rightmost column of Table I we 
exhibit the values of F 2 /£, with all the numerals in the lattice unit. 

Note that in the parameter region 
VI7|m| 3> 1, eq.Q should approach the 
a model of nonchiral parametric RM en- 
sembles Jl3 but the pion decay or diffusion 
constant multiplying tr BWBU is unaf- 
fected. Accordingly, if the mean level 
spacing is approximately constant in a 
window in the very vicinity of the ori- 
gin, one can determine F 2 / E from the 
bulk correlation (namely the LSD) in that 
window. As LSDs admit a high precision 
determination of p and achieve extremely 
high linearity of the ip-y/Ap plots than 




Fig. 4. Low-energy constants F 2 and E for 
SU(2) quenched lattice gauge theory. 



those from SEDs (Fig. 3), we have adopted the former for the determination of F 2 / E. 
Coupling dependence of these low-energy constants are summarized in Table I and 
plotted in Fig. 4 (including the result from simulations at (3 = 1.5). In order to offset 
the number of components in the gauge multiplet, these constants are multiplied by 
2/3 for adjoint in the plot. At (3 = 0, both low-energy constants agree between fun- 
damental and adjoint representations. This observation is consistent with the fact 
that in the strong coupling limit these constants are common to the large- A Sp(2A) 
and 0(A) lattice gauge theories JIB 1 which share the same antiunitary symmetries as 
SU(2)-fundamental and adjoint. 



Summary. We have evaluated the smallest eigenvalue and level spacing distribu- 
tions for (ch)GSE-(ch)GUE and (ch)GOE-(ch)GUE crossover using the Nystrom- 
type method, the former being our new contribution. These RM results are applied 
to fit the fundamental and adjoint staggered Dirac spectra of SU(2) quenched lattice 
gauge theory under twisted boundary condition. Excellent one-parameter fitting is 
achieved for all cases of concern. Acute sensitivity of our fitting distributions, pi(s) 
and P(s), on the crossover parameter p leads to the precise determination of the 
pion decay constant F from its twisting dependence. This method, feasible on a 
small-size lattice, has a clear advantage over the conventional method using axial 
correlators, which inevitably requires a large temporal dimension. 

Our treatment is complementary to the previous approach of determining F of 
two-color QCD from its Dirac spectrumJM that conversely measured the response 
of complex Dirac eigenvalues to a real chemical potential. In a practical point of 
view, the use of imaginary chemical potential has an advantage because (i) it needs 
no projection of eigenvalues to the real or imaginary axis for fitting, which is usually 
requisite for dealing with the complex eigenvalues, and (ii) two-dimensional motion of 
complex eigenvalues may lead to large statistical fluctuation. Similar complementary 
treatments were applied for three-color QCD at real and imaginary isospin chemi- 
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cal potential, corresponding to non-Hermitian chiral and Hermitian crossover 

chiral RMsp^'H^} respectively. Combined with the results reported here that fills the 
vacancy, the established fact that Dirac spectra in all three cases (SU(2)-fund.+ / u, 
SU(2)-adj.+/U, SU(3)-fund.+^i so ) agree perfectly with predictions from correspond- 
ing zero-mode-approximated chiral Lagrangians in both regions of fi 2 ^ constitutes 
an encouraging evidence for the validity of analytic continuation in the //-plane. 

Finally we note that our preliminary simulation have affirmed that Dirac spectra 
of SU(2)xU(l) quenched lattice gauge theory also fit excellently to the parametric 
RM predictions. We are currently accumulating numerical data on lattices of larger 
size than treated here, and the results will be reported in a subsequent publication. 
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